Statistically Preserved Structures and Anomalous Scaling in Turbulent Active Scalar Advection

The anomalous scaling of correlation functions in the turbulent statistics of active scalars (like temperature in turbulent convection) is understood in terms of an auxiliary passive scalar which is advected by the same turbulent velocity field. While the odd-order correlation functions of the active and passive fields differ, we propose that the even-order correlation functions are the same to leading order (up to a trivial multiplicative factor). The leading correlation functions are statistically preserved structures of the passive scalar decaying problem, and therefore universality of the scaling exponents of the even-order correlations of the active scalar is demonstrated.Comment: 4 pages, 5 figures, submitted to Phys. Rev. Let

Self-organising management of Grid environments

This paper presents basic concepts, architectural principles and algorithms for efficient resource and security management in cluster computing environments and the Grid. The work presented in this paper is funded by BTExacT and the EPSRC project SO-GRM (GR/S21939)

Cumulus cloud venting of mixed layer ozone

Observations are presented which substantiate the hypothesis that significant vertical exchange of ozone and aerosols occurs between the mixed layer and the free troposphere during cumulus cloud convective activity. The experiments utilized the airborne Ultra-Violet Differential Absorption Lidar (UV-DIAL) system. This system provides simultaneous range resolved ozone concentration and aerosol backscatter profiles with high spatial resolution. Evening transects were obtained in the downwind area where the air mass had been advected. Space-height analyses for the evening flight show the cloud debris as patterns of ozone typically in excess of the ambient free tropospheric background. This ozone excess was approximately the value of the concentration difference between the mixed layer and free troposphere determined from independent vertical soundings made by another aircraft in the afternoon

Probability Density Function of Longitudinal Velocity Increment in Homogeneous Turbulence

Two conditional averages for the longitudinal velocity increment u_r of the simulated turbulence are calculated: h(u_r) is the average of the increment of the longitudinal Laplacian velocity field with u_r fixed, while g(u_r) is the corresponding one of the square of the difference of the gradient of the velocity field. Based on the physical argument, we suggest the formulae for h and g, which are quite satisfactorily fitted to the 512^3 DNS data. The predicted PDF is characterized as (1) the Gaussian distribution for the small amplitudes, (2) the exponential distribution for the large ones, and (3) a prefactor before the exponential function for the intermediate ones.Comment: 4 pages, 4 figures, using RevTeX3.

Optimized Broadcast for Deep Learning Workloads on Dense-GPU InfiniBand Clusters: MPI or NCCL?

Dense Multi-GPU systems have recently gained a lot of attention in the HPC arena. Traditionally, MPI runtimes have been primarily designed for clusters with a large number of nodes. However, with the advent of MPI+CUDA applications and CUDA-Aware MPI runtimes like MVAPICH2 and OpenMPI, it has become important to address efficient communication schemes for such dense Multi-GPU nodes. This coupled with new application workloads brought forward by Deep Learning frameworks like Caffe and Microsoft CNTK pose additional design constraints due to very large message communication of GPU buffers during the training phase. In this context, special-purpose libraries like NVIDIA NCCL have been proposed for GPU-based collective communication on dense GPU systems. In this paper, we propose a pipelined chain (ring) design for the MPI_Bcast collective operation along with an enhanced collective tuning framework in MVAPICH2-GDR that enables efficient intra-/inter-node multi-GPU communication. We present an in-depth performance landscape for the proposed MPI_Bcast schemes along with a comparative analysis of NVIDIA NCCL Broadcast and NCCL-based MPI_Bcast. The proposed designs for MVAPICH2-GDR enable up to 14X and 16.6X improvement, compared to NCCL-based solutions, for intra- and inter-node broadcast latency, respectively. In addition, the proposed designs provide up to 7% improvement over NCCL-based solutions for data parallel training of the VGG network on 128 GPUs using Microsoft CNTK.Comment: 8 pages, 3 figure

Real-time Loss Estimation for Instrumented Buildings

Motivation. A growing number of buildings have been instrumented to measure and record earthquake motions and to transmit these records to seismic-network data centers to be archived and disseminated for research purposes. At the same time, sensors are growing smaller, less expensive to install, and capable of sensing and transmitting other environmental parameters in addition to acceleration. Finally, recently developed performance-based earthquake engineering methodologies employ structural-response information to estimate probabilistic repair costs, repair durations, and other metrics of seismic performance. The opportunity presents itself therefore to combine these developments into the capability to estimate automatically in near-real-time the probabilistic seismic performance of an instrumented building, shortly after the cessation of strong motion. We refer to this opportunity as (near-) real-time loss estimation (RTLE). Methodology. This report presents a methodology for RTLE for instrumented buildings. Seismic performance is to be measured in terms of probabilistic repair cost, precise location of likely physical damage, operability, and life-safety. The methodology uses the instrument recordings and a Bayesian state-estimation algorithm called a particle filter to estimate the probabilistic structural response of the system, in terms of member forces and deformations. The structural response estimate is then used as input to component fragility functions to estimate the probabilistic damage state of structural and nonstructural components. The probabilistic damage state can be used to direct structural engineers to likely locations of physical damage, even if they are concealed behind architectural finishes. The damage state is used with construction cost-estimation principles to estimate probabilistic repair cost. It is also used as input to a quantified, fuzzy-set version of the FEMA-356 performance-level descriptions to estimate probabilistic safety and operability levels. CUREE demonstration building. The procedure for estimating damage locations, repair costs, and post-earthquake safety and operability is illustrated in parallel demonstrations by CUREE and Kajima research teams. The CUREE demonstration is performed using a real 1960s-era, 7-story, nonductile reinforced-concrete moment-frame building located in Van Nuys, California. The building is instrumented with 16 channels at five levels: ground level, floors 2, 3, 6, and the roof. We used the records obtained after the 1994 Northridge earthquake to hindcast performance in that earthquake. The building is analyzed in its condition prior to the 1994 Northridge Earthquake. It is found that, while hindcasting of the overall system performance level was excellent, prediction of detailed damage locations was poor, implying that either actual conditions differed substantially from those shown on the structural drawings, or inappropriate fragility functions were employed, or both. We also found that Bayesian updating of the structural model using observed structural response above the base of the building adds little information to the performance prediction. The reason is probably that Real-Time Loss Estimation for Instrumented Buildings ii structural uncertainties have only secondary effect on performance uncertainty, compared with the uncertainty in assembly damageability as quantified by their fragility functions. The implication is that real-time loss estimation is not sensitive to structural uncertainties (saving costly multiple simulations of structural response), and that real-time loss estimation does not benefit significantly from installing measuring instruments other than those at the base of the building. Kajima demonstration building. The Kajima demonstration is performed using a real 1960s-era office building in Kobe, Japan. The building, a 7-story reinforced-concrete shearwall building, was not instrumented in the 1995 Kobe earthquake, so instrument recordings are simulated. The building is analyzed in its condition prior to the earthquake. It is found that, while hindcasting of the overall repair cost was excellent, prediction of detailed damage locations was poor, again implying either that as-built conditions differ substantially from those shown on structural drawings, or that inappropriate fragility functions were used, or both. We find that the parameters of the detailed particle filter needed significant tuning, which would be impractical in actual application. Work is needed to prescribe values of these parameters in general. Opportunities for implementation and further research. Because much of the cost of applying this RTLE algorithm results from the cost of instrumentation and the effort of setting up a structural model, the readiest application would be to instrumented buildings whose structural models are already available, and to apply the methodology to important facilities. It would be useful to study under what conditions RTLE would be economically justified. Two other interesting possibilities for further study are (1) to update performance using readily observable damage; and (2) to quantify the value of information for expensive inspections, e.g., if one inspects a connection with a modeled 50% failure probability and finds that the connect is undamaged, is it necessary to examine one with 10% failure probability

Turbulent Drag Reduction by Flexible and Rodlike Polymers: Crossover Effects at Small Concentrations

Drag reduction by polymers is bounded between two universal asymptotes, the von-K\'arm\'an log-law of the law and the Maximum Drag Reduction (MDR) asymptote. It is theoretically understood why the MDR asymptote is universal, independent of whether the polymers are flexible or rodlike. The cross-over behavior from the Newtonian von-K\'arm\'an log-law to the MDR is however not universal, showing different characteristics for flexible and rodlike polymers. In this paper we provide a theory for this cross-over phenomenology.Comment: 5 pages, 4 figures, submitted to Physical Review

Dynamics of Scalar Fields in the Background of Rotating Black Holes

A numerical study of the evolution of a massless scalar field in the background of rotating black holes is presented. First, solutions to the wave equation are obtained for slowly rotating black holes. In this approximation, the background geometry is treated as a perturbed Schwarzschild spacetime with the angular momentum per unit mass playing the role of a perturbative parameter. To first order in the angular momentum of the black hole, the scalar wave equation yields two coupled one-dimensional evolution equations for a function representing the scalar field in the Schwarzschild background and a second field that accounts for the rotation. Solutions to the wave equation are also obtained for rapidly rotating black holes. In this case, the wave equation does not admit complete separation of variables and yields a two-dimensional evolution equation. The study shows that, for rotating black holes, the late time dynamics of a massless scalar field exhibit the same power-law behavior as in the case of a Schwarzschild background independently of the angular momentum of the black hole.Comment: 14 pages, RevTex, 6 Figure

Unconventional Gravitational Excitation of a Schwarzschild Black Hole

Besides the well-known quasinormal modes, the gravitational spectrum of a Schwarzschild black hole also has a continuum part on the negative imaginary frequency axis. The latter is studied numerically for quadrupole waves. The results show unexpected striking behavior near the algebraically special frequency $\Omega=-4i$. This reveals a pair of unconventional damped modes very near $\Omega$, confirmed analytically.Comment: REVTeX4, 4pp, 6 EPS figure files. N.B.: "Alec" is my first, and "Maassen van den Brink" my family name. v2: better pole placement in Fig. 1. v3: fixed Refs. [9,20]. v4: added context on "area quantum" research; trimmed one Fig.; textual clarification

Wave Propagation in Gravitational Systems: Completeness of Quasinormal Modes

The dynamics of relativistic stars and black holes are often studied in terms of the quasinormal modes (QNM's) of the Klein-Gordon (KG) equation with different effective potentials $V(x)$. In this paper we present a systematic study of the relation between the structure of the QNM's of the KG equation and the form of $V(x)$. In particular, we determine the requirements on $V(x)$ in order for the QNM's to form complete sets, and discuss in what sense they form complete sets. Among other implications, this study opens up the possibility of using QNM expansions to analyse the behavior of waves in relativistic systems, even for systems whose QNM's do {\it not} form a complete set. For such systems, we show that a complete set of QNM's can often be obtained by introducing an infinitesimal change in the effective potential